made only after the accumulation of such a large 
amount of data as this. The small peaks at the 
right of figure 3 indicate a slight tendency for the 
observed counts to be influenced by effects other 
than statistical fluctuations. However, it must be 
remembered that the data include instances where 
either counters or scaling circuits have gone bad 
and have evidenced a need for repair or replacement, 
so that this slight excess is to be expected. It is 
believed that this is the first time that such a large 
collection of Chi square values has been obtained 
from the ratio of the error variance of Latin squares 
to the expected (Poisson) variance. 

It is well known that individual measurements of 
radioactive disintegrations follow a Poisson law, 
but the close agreement between the experimental 



and theoretical curves illustrates one further point: 
that the Latin square arrangement has been com- 
pletely successful in eliminating extraneous, non- 
Poisson fluctuations from the measurements. 



The author thanks W. J. Youden and J. M. 
Cameron, of the Statistical Engineering Section, for 
their suggestions as to the use of the Latin squares 
and for many valuable hours of discussion; and 
Margaret Selgin and Lucy Cavallo, of the Radio- 
activity Section, who performed the calculations 
summarized in figure 3. 

Washington, July 26, 1950. 
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A Problem in Precision Cam Design 

By Joseph Blum 

This paper proposes an analytic method for determining the profile of a cam required in 
a device that demands high precision. The method is applicable whenever the equation for 
the profile can be obtained by making a less stringent assumption, namely, that the follower 
makes contact at a fixed point of its extremity. The true profile can then be determined 
numerically by considering it as the envelope of a family of curves. 



The design of an X-ray spectrograph often requires 
that a constant rate of oscillation be maintained in 
the crystal holder. The crystal rotates through an 
angle # (the angle of travel), reverses its direction 
and then proceeds with the same angular velocity. 
This type of motion may be obtained from a cam 
rotating with constant angular velocity, co. An 
analytic method for computing the profile of the 
cam is presented in this paper which extends, in a 
certain sense, the application of a paper written by 
J. B. Friauf; 1 the notation therein will be preserved 
to a considerable extent. 



i J. B. Friauf, The design of a cam for an X-ray spectrograph, J. Opt. Soc. Am" 
Rev. Sci. Instr. 11, 289 to 296 (1925). 




In figure 1, a, represents the axis about which 
the cam rotates, A the axis about which the crystal 
holder rotates, and B the center of the ball at the 
end of the follower. The follower ACB makes 
contact with the cam and serves to rotate the crystal 
holder. The distances AO and AB are equal; denote 
this common length by B. Let T\ and r 2 be, respec- 
tively, the least and greatest radii measured from 
to the point B; let r be the variable radius, OB 
corresponding to the angle a, which r makes with the 
least radius r\. 

Denote the time by t and let £=0 represent the 
time when the least radius coincides with the line 
segment extending from to A. Then the angle 




Figure 1. 
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between r x and OA is given by ut. Let denote 
angle OAB and </> the angle AOB. Unless otherwise 
indicated, all angles are in radians. 
From figure 1, a, it is clear that 



sin 



r 

~2R' 



(1) 



(2) 



from which it follows that the smallest angle 0i 
occurs when r=r x and the largest angle 2 when 
r=r 2 , or 

* m 2=2R> 
and 

02 ^2 /r>\ 

Sm 2=2S (3) 

It is also clear that $ = B 2 — B x . Differentiating eq 1 
yields 

dd 1 dr 1 dr da 



C °*2dt Rdt Rdadt 

bince a=coc— (/>, -jt— 

Substituting this into eq 4 yields 
6 dd 1 dr 



(4) 



d$ 
' dt 



cos 



2 dt R d 



r ( dcf>\ 
ar~Tt)' 



(5) 



Again, since 0=-(tt— 0), (d(t>/dt= 
stituting this into eq 5 results in 



■%(d$/dt) 9 and sub- 



cos 



d0_J_ dr 
2 d* ~7* 



dr / 1 dd\ 
daV + 2dtf 



(6) 



The follower rotates through the angle /3 during the 
same time interval that it takes the cam to rotate 
through t radians. Therefore, as the respective 
angular velocities are constant, 



d0 

dt d6 0co 

CO 7T at TT 



Substituting this into eq 6 gives 



cos 



>_:[ dr Ar 1\ 
i~~Rda\P + 2j' 



From eql cos (d[2) = (l/2R)-y/±R 2 -r 2 , 
this into eq 7 and simplifying yields 



(7) 
Substituting 



2 dr 



dc 



V4# 2 



/Si 



(8) 



where j8i = tt//3 + & 

After integrating this equation there is obtained 

where C\ is the constant of integration. 



(9) 



When a=0, r=r x ) solving for C\ results in Ci = 2 
sin ~ 1 ri/2R=di. Equation 9 can be expressed in the 
form 



r=2R sin - 



2(j8 



;+0| 



(10) 



The greatest radius, r 2 , occurs when a = Tr J rP/2, 
which may be vertified by substituting into eq 10. 

If the ball that makes contact with the cam were 
sufficiently small, then eq 10 would give one part of 
the cam's contour, and a similar equation would give 
the remaining part. Frequently, however, this con- 
dition is not met, and the method given in the ref- 
erence in footnote 1 must be extended to secure a 
high level of precision. Consider then the family of 
circles of radius p, whose centers lie on the curve 
given by eq 10. The envelope of this family con- 
sists of two distinct curves; the one that lies nearer 
the center of curvature of eq 10 coincides with the 
contour of the cam. This envelope can be obtained 
by a numerical determination of the ^-discriminant. 2 

Let the curve of eq 10 be given in Cartesian co- 
ordinates by 



y=f(x). 



(ii) 



Let (x,y) be an arbitrary point on this curve. The 
equation of the circle with radius p and center (x,y) 
is given by 

(X-xy+(Y-yY=p\ (12) 

Substituting y=j(x) into eq 12, there is obtained a 
one-parameter family of circles 

{X-xY+{Y-j{x)Y= P \ (13) 

Differentiating eq 13 with respect to x yields 

{X-x}+f(x)lY-f(x)}=:0. (14) 

The ^-discriminant is obtained by solving eq 13 and 
14 simultaneously for a discrete set of values of the 
parameter x. The method outlined below yields the 
numerical solution: 

1. For a selected set of values a t ) determine the 
values of r t by using eq 10. Denote these points 
by P.. 

2. Let n be the angle of inclination of the tangent 
line at P*; let ^ be the angle between the radius r t 
and the tangent line. Then, since r i =a i J r\(/ i , 



/ / (x) = tanr i = 



tanc^ + tan^ 
1 — tan on tan \f/ t 



Now tan ypi=r i jr , i) where r/ is the derivative dr/da 
evaluated at P t . This leads to 

tan^ = 2/3 1 tani(^+0 1 ). 



■2\0« 



2 E. L. Ince, Ordinary differential equations, p. 85 to 86 (Longmans, Green 
& Co., New York, N. Y., 1944). 
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These relations enable one to compute the value of 

3. The solution of eq 13 and 14 is easily obtained 
in the form 



Yi= yi 



Vi + i/'fe)} 2 
p 



(15) 



Vi + i/'fe)} 2 

The points (X t , Y t ) are on the envelope, as required. 

To obtain the other part of the cam, consider that 
phase of the motion in which the angle 6 decreases 
(fig. 1, b). For this phase, 



dd 
dt 



0a> 

IT 



Substituting this into eq 6 results in 
6 1 dr /l tt\ 



(16) 



where 7 denotes the angle that r makes with the 
greatest r adius r 2 (see fig l,b). Substituting cos (6/2) = 
(1/2R) ^4:R 2 —r 2 into eq (16) and simplifying gives 



2dr 



dy 



V4^ 2 -r 2 £ 2 
where /3 2 =l/2 — ir/0. 

When eq 17 is integrated there results 



2sin " 1 i=i+^ 



(17) 



(18) 



where C 2 is the constant of integration. When 7=0, 
r=r 2 ; solving for C 2 yields C 2 =2 sin _1 r 2 /2i?=0 2 . 
Equation 18 can be expressed in the form 



r=2R sin 



3te+*> 



(19) 



As a check, when y = w— /3/2, r=2R sin % (d 2 —p) = 




Figure 2. 

circles with centers on the curve (19) and radius p 
leads to the remainder of the contour of the cam. 

To complete the analysis, consider figure 2. The 
arc ABC is the locus of eq 10, and arc A'RB'C is 
its envelope. Similarly, the arc DEF is the locus of 
eq 19; arc D'RE'F' is its envelope and may be found 
numerically by following the procedure outlined 
above for eq 10. The high point of the cam is 
located numerically by solving for the intersection 
of the two envelopes. The low portion of the cam 
does not coincide with the envelope, as is evident 
from the figure. Let T be the point of intersection 
of the curves 10 and 19. The circular arc, PQ, with 
center at T and radius p constitutes the low portion 
of the cam. These considerations enable one to 
complete the numerical solution of the cam design. 

In summary, it should be noted that the essential 
idea in this approach is to first find the equation of 
the contour, assuming that the follower makes con- 
tact at a fixed point of its extremity; then, by using 
the true form of the follower, obtain the actual cam 
contour by determining the envelope of a family of 
curves. It is evident that the analytic method ap- 
plied in this particular case may be used in other 
problems of precise cam design. 



2R sin y 2 



6 1 =r 1 . The envelope of the family of Washington, August 4, 1950. 
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